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Abstract 

In the context of AdS/CFT with flavor, we consider the type IIB supergravity solution cor- 
responding to a fully localized D3/D7 intersection. We complete the standard metric ansatz by 
providing an analytic expression for the warp factor, under the assumption of a logarithmically 
running axion-dilaton. From the gauge dual perspective, this behavior is related to the positive 
beta function of the N = 4, d = 4 SU (N c ) super Yang- Mills gauge theory, coupled to Nf funda- 
mental N = 2 hypermultiplets. We comment on the existence of tadpoles and relate them to the 
same gauge theory beta function. Next we consider a classical spinning string configuration in 
the decoupling limit of the D3/D7 geometry and extract the flavor (Nf) dependence of the asso- 
ciated meson Regge trajectory. Including the backreaction of the D7-branes in the supergravity 
dual allows for going beyond the quenched approximation on the dual gauge theory side. 



1 Introduction 



In the recent past, much effort has gone into studying the AdS/CFT duality of Yang-Mills 
theories with matter in the fundamental representation of the gauge group. Starting with 
the work of Karch and Katz pQ, fundamental matter was introduced as a new open 
string sector arising through the embedding of probe branes into various supergravity 
backgrounds [2 -[24j. This approach corresponds to the quenched approximation in lattice 
QCD which ignores the effect of the creation and annihilation of virtual quark- ant iquark 
pairs. 

In a parallel development, supergravity backgrounds of fully localized Dp/D(p+4) 
brane intersections have been constructed which include the backreaction of the "flavor" 
D(p+4) brane [213 [2E1 123 12E1 [2312111 EH]- Such brane configurations allow for the discussion 
of flavors beyond the quenched approximation. An interesting example is the D3/D7 
brane intersection which, in contrast to the same system in the probe approximation [T|, 
includes the backreaction of the D7-branes on the geometry. This configuration describes 
an Af = 2 gauge theory obtained from the coupling of M = 4 SU(N C ) super Yang- 
Mills theory to Nf hypermultiplets in the fundamental representation of the gauge group. 
Similar to QED and 4 theory, the theory has a positive beta function leading to an 
ultraviolet divergence in the gauge coupling constant. An ultraviolet Landau pole usually 
indicates the breakdown of the perturbative field theory at high energies, and it is not yet 
understood how this pathology can be cured. 

In this paper we complete the construction of the type IIB supergravity solution for 
the D3/D7 intersection restricted to the near-core region of the D7-branes, where the 
axion-dilaton has a simple logarithmic behavior. Unlike the full complex structure of the 
D7-branes, which is finite everywhere, the dilaton runs into a singularity at some finite 
distance away from the D7-branes. 

As we will show, there is a one-to-one map between the logarithmic dilaton profile 
and the one-loop running of the gauge coupling. This precisely maps the Landau pole 
in the field theory to the dilaton divergence in the supergravity solution. Similarly, the 
chiral anomaly of the U(1)k symmetry, which leads to a non-vanishing theta angle, is 
appropriately represented by a non-trivial axion field. 

As shown in |28| l29j. the construction of a supergravity solution for the D3/D7 system, 
reduces to the problem of finding a solution to a Poisson equation for the warp factor. 
Despite claims to the contrary, it is possible to find an analytical solution for the warp 
factor. We show that the Poisson equation can be reduced to a linear ordinary differential 
equation which can be solved by means of techniques developed by Gesztesy and Pittner 
|32j . The general solution is represented by a uniformly convergent series expansion for 
which we give an explicit proof of convergence. 

Due to the logarithmic ultraviolet divergence, the string background is not finite and 
is expected to have an uncanceled tadpole. Nevertheless, as pointed out in [33], if the 
tadpole is logarithmic, the background is still consistent as far as the embedding of the 
non-conformal gauge theory is concerned. Logarithmic tadpoles do not reflect gauge 
anomalies, but instead provide the correct one-loop running of the gauge coupling. Using 
results found in [33], we will show that the open string one- loop annulus amplitude has 
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indeed the correct logarithmic behavior expected from the positive beta function of the 
M = 2 field theory. 

The universal structure of the D3/D7 configuration gives rise to many applications 
in string theory. Apart from its application to strong interactions, the D3/D7 system 
has also received renewed attention in a cosmo logical context. A prominent example is 
the D3/D7 model in [33] which has an effective description as hybrid inflation. It might 
even turn out that D3 and D7-branes are the preferred branes in brane cosmology, as can 
be argued by a discussion of the decay and annihilation rates of D-branes in a higher- 
dimensional universe [SH]. Finally, the D3/D7 plays an essential role in the construction 
of semi- realistic four-dimensional string vacua [3*7] . 

In the second part of the paper, we study meson Regge trajectories in the Af = 2 
gauge theory dual to the previously constructed D3/D7 supergravity solution. We are 
particularly interested in the dependence of the meson spectrum on the number of sea 
quark flavors introduced by the D7-branes. 

Following the analysis in [3 El OH] > we determine the spectrum of mesons with large 
spin J by computing the energy and angular momentum of a semi-classical string rotating 
in the near-horizon geometry of the D3/D7 intersection. This string is attached to an 
additional probe D7-brane which corresponds to a massive flavor Q in the background of 
Nf massless flavors qi. The setup is similar to the one in but with Nf flavors of light 
sea quarks turned on. 

By varying the number of D7-branes Nf and keeping the mass of the single flavor Q 
fixed, we determine the Regge trajectory of the meson QQ. We find that for a meson with 
a particular spin J, the energy as well as the string tension decreases with the number 
of massless flavors, whereas the string length, which measures the distance between the 
quarks inside the meson, remains unaffected by the sea quark flavors. 

In the field theory, the decrease in the string tension can be ascribed to the influence 
of virtual color-singlet qq pairs on the color flux inside the meson QQ. As it is well- 
known, virtual qq pairs cause a screening of the color charge, effectively diminishing 
the color force and with it the string tension. This quark-loop effect amounts to going 
to first order in Nf/N c beyond the quenched approximation in the gauge theory and 
arises as a consequence of including the backreaction of the D7-branes in the D3/D7 
supergravity background. According to the standard lore of AdS/CFT, a gauge theory 
one-loop quantum effect is captured at the classical level by the supergravity dual. 

The paper is organized as follows. In Sec. 121 we revisit the D3/D7 configuration and 
derive an analytical expression for the warp factor in terms of a uniformly convergent se- 
ries expansion. From this we also deduce the supergravity solution for fractional D3/D7 
branes which describes M = 2 SU(N C ) super Yang-Mills theory with Nf fundamental 
hypermultiplets, but without an adjoint hypermultiplet. We further discuss several issues 
related to the consistency of the string theory embedding of a non-conformal gauge theory 
with running gauge coupling. In Sec. El we discuss the flavor dependence of Regge trajec- 
tories in the D3/D7 brane theory. A summary of the D3/D7 supergravity solution can be 
found in the conclusions in Sec. EJ App. [X] contains a review on the D7-branes geometry. 
App. [B] contains the convergence proof of the series expansion of the warp factor. 
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2 The fully backreacted D3/D7 solution 



In the following we construct the supergravity solution for a fully localized D3/D7 in- 
tersection in flat space. The corresponding near-horizon limit of this solution is dual to 
M = 2 super Yang-Mills theory with Nf hypermultiplets in the fundamental and one 
hypermultiplet in the adjoint representation of the gauge group SU(N C ). This set-up is 
T-dual to the D2/D6 and D4/D8 geometry derived in j2S] and [2B|, respectively. Previous 
work on the D3/D7 solution in flat space can be found in [23 12H1 129*1 15D] . 

2.1 The D3/D7 solution 

The D3/D7 brane intersection consists of a stack of iV c coincident D3-branes which is 
embedded into the world volume of Nf D7-branes. The D7-branes are located in ten- 
dimensional flat space such that they share four longitudinal directions with the D3-branes 
while wrapping or spanning four out of the six transverse directions. This embedding may 
be represented pictorially as 
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and preserves 8 supersymmetries as well as a SO (A) x 5*0(2) isometry. Note that sepa- 
rating the D3-branes from the D7-branes in the 89 direction would explicitly break the 
S0(2). 

Following [2H1 HUE ED] > we make the general metric ansatz 

ds 2 w = hT x l\y m )dxl + h l l 2 (y m )g mn dy m dy n , (2.1) 

where x M (/i = 0, . . . , 3) are coordinates on the longitudinal spacetime, and y m (m = 
4, . . . , 9) are coordinates transverse to the D3-branes. 

The six-dimensional Kahler metric transverse to the D3-branes is completely fixed by 
the D7-brane metric 

ds 2 6 = dzidzi + dz 2 dz 2 + e 9 ^'^ dz 3 dz 3 , (2.2) 
where the function e* is given by 

=r 2 (z 3 )\ V (r)\ i \z 3 \~ N f/ 6 . (2.3) 

For a review on the properties of the D7-brane metric, see App. ["""] 

We now have to make a choice for the complex structure r. The full D3/D7 solution 
would require the (full) complex structure r of the D7-branes as given by Eq. (O in 
the appendix. However, for reasons which become apparent in Sec. 12.51 we consider 
only the weak coupling region of the D7-branes, in which the complex structure is well- 
approximated by 

r \ ■ N f^ ( Zs 

27T \p L 
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The integration constant pL is given by 



■>TT 



PL = e^ N f , (2.5) 

such that = g s for Nf = 0. 

This corresponds to focusing in the region close to the D7-branes, \z\ <C Pl, where e* 
simplifies to 

e *(*a.*0 „ T2 ( Z3 ) = I _ Hi logZ3 ^ (2 .6) 

It has been shown in j2H] that this ansatz preserves the right amount of supersymmetry. 
What remains to be done is to derive an analytical expression for the warp factor h(y m ) 
in the metric ansatz (|2.1|) . 

2.2 The warp factor 

In order to find an analytical expression for the warp factor, we must solve the transverse 
Laplacean 

n (6) h=(d 1 d 1 + d 2 d s + e- 9 d 3 d s )h. (2.7) 

This expression for the warp factor was first examined in j2H], where it was solved to a 
first order approximation z 3 = Z 3 + 5z 3 around an arbitrary but fixed point Z 3 away from 
the D7-branes. 

We next introduce real coordinates on the space transverse to the D3-branes, 

r 2 = \z 1 \ 2 + \z 2 \ 2 , z 3 = pe^, (2.8) 
in terms of which the transverse metric ()2.2|) becomes 

ds 2 6 = dfdr + e*(dp 2 + p 2 dip 2 ) . (2.9) 
In these coordinates the Laplacean reads 

D w h(r, p, if) = ^d r r% + e~ 9 ^ (~d p pd p + L&\ ^ h (r, p, <p) . (2.10) 

The warp factor is, in fact, the Green's function of the Laplace equation defined by 

D 6 G(p, ^ r; p', if', P) = N c e~*5 4 (r- P)5{p - p')-5(tp - if') . (2.11) 

P 

We solve for the Green's function by performing an expansion in Fourier modes 

p) 5 (cp - V ') = -L-J d'qe^^ ( 2 - 12 ) 
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and similarly 

G(p, <p, f; p\ <p', r>) = 1 + Q m I d'qe^-^M pV W) , (2-13) 

l 

where Qdz = 4ng s NJs is the D3-brane charge. Substituting both the expansion of the 
Green's function and that of the delta function into (j2.11j) . we find 

- ^e-V(p;p') + e-"-d p pd pm , g (p;p') - q 2 m, g (p; p) = j^^-Ap- p') ( 2 - 14 ) 
p p i z7r j p 

with q = \q\. 

Then for p ^ 0, Eq. ()2.14j) leads to the differential equation 

d 2 1 d l 2 \ 

+ V(p) + -) yi , q (p;p')=0, (2.15) 



dp 2 p dp p 2 
with the logarithmic potential 



1 Nf^Aj* 



V(p)= ---^logp q 2 . (2.16) 
\9s 2vr J 

By means of the redefinition of the radial coordinate 

x = log p/ p L , (2.17) 
the differential equation (|2.15|) turns into 

(d 2 x - l) yi , q {x) = \xe 2x yUx) , A = ^pW , (2.18) 

where we defined y(x) = y(p(x);p' = 0). We note that in the new coordinate x, the 
near-core region (p Pl) is located near to x — > — oo. 

There are two independent solutions to the differential equation (|2.18|) : they are dis- 
tinguished by their behavior near the origin and at infinity. Let us denote by yi, q (p) the 
solution that is well behaved near the origin (if Nj = 0, this solution reduces to the mod- 
ified Bessel function Ii(qp)), and by yi, q (p) the solution which is well behaved at infinity 
(similarly, for Nf = 0, this solution becomes the modified Bessel function Ki(qp)). The 
explicit construction of these two independent solutions will be addressed at length in the 
paragraphs below. The Green's function G(p, ...;//,...) is then obtained by taking the 
product yi, q (p<)yi,q(p>), where p< = min(p, p') and p> = max(p,p'). 

Without loss of generality we proceed to construct the Green's function G(r, p, <p; 0) 

G(r,p,<p;0) = l + Q m r^~ 3 q 3 f dty/T=Pe«** £ y l>q (0) yi Me il * 

^° 71 l=-oo 

= 1 + Qd3 J" dq iqr) ^ r) y , q (p) , (2.19) 

in which case we use that 2/o,g(0) is the only non- vanishing term. This restricts the sum 
over I to one term corresponding to I = 0. 
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Near-core solution of Eq. (I2.18R 

Before deriving the general solution of ([2. 18)1 . let us first consider the solution in the 
near-core region of the D7-branes. Here it is convenient to define the coordinate 

p 2 = \xe 2x . (2.20) 

The corresponding differentials satisfy 

2pdp = Xe 2x (l + 2x)dx. (2.21) 

In the near-core region, close to x — > — oo, we may neglect the term \e 2x dx on the 
right hand side of (J2.18j) and obtain 

dp^pdx. (2.22) 
In this approximation, Eq. (|2.18|) reduces to the differential equation (y = y 0tq ) 

fatfM=iWD, (2-23) 

whose general solution is 

y{p) = Cl I (p) + C2K Q (p) , (2.24) 

where C\ and c<i are constants, and Iq and K are the modified Bessel functions of the first 
and second kind, respectively. Since we are interested in solutions which are singular at 
the location p = of the D7-branes and well behaved at infinity, we set C\ = 0, C2 = 2n 2 
and obtain 

y(p) = 2tt 2 K (p) . (2.25) 
Substituting this solution into (j2.19j) . we find the near-core harmonic function 



h(r,p) = 1 + Q D3 !°° dqq 2 ^^-K Q {V\x~^) 
Jo 2r 



" (r^ - p^ log JL)* + (r» + p 2 e*f ' ^ 

V r 2lT PL 

where we used the identity 

\ xe 2x = p 2 q 2 (L-^l\ gp]= p 2 g 2 e* . (2.27) 
\9a 2vr ) 

For Nf = 0, we recover the harmonic function of D3-branes as required. 

The near-core solution can be generalized to the case when the D7-branes are separated 
from the D3-branes by a distance d in the overall transverse plane. Then, the function 
e* has the form 

e*^) = 1 _ ^ log(p 2 - 2pd cos <p + d 2 ) . (2.28) 
The distance d corresponds to the mass m of the Nf flavors. 
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Exact solution of Eq. (IQl 

Studying the Schrodinger equation of electrons in a logarithmic potential, Gesztesy and 
Pittner (GP) [32] found a particular solution to (j2.18|) which is however not the most 
general one. As shown in Appendix El the GP solution is not of immediate use for the 
D3/D7 system, since it reduces to Io(pq) rather than K (pq) in the absence of D7-branes. 
It is therefore not possible to recover the harmonic function of D3-branes in the limit 
Nf = 0. In the following, we construct a solution for a finite number of flavors, Nf ^ 0, 
which reduces to the D3-brane solution for Nf = 0. 

Let us first consider the asymptotic behavior of solutions to ()2.18|) . At x — > — oo, there 
are two different solutions, 

y(x) — c and y(x) = ax + b , (2.29) 

with arbitrary constants a, b, c (a 7^ 0). The GP solution [32| behaves like y(x) = Iq{0) = 1 
at x — > — 00. We are however interested in solutions with a singular behavior at x — > —00. 
In particular, the sought after solution should behave as 

y{p) = K (pq) w - log(pg/2) - 7 . (2.30) 

As discussed in Appendix [BJ this is equivalent to choosing the boundary condition y(x) = 
ax + b with 

a = — 1 , b = — x — 7 , (2.31) 

where Xq = logpig/2 and 7 = —^(1) is the Euler-Mascheroni constant. 

Following [321, the solutions with this asymptotic behavior are given by the series 

00 

y(x) = 27T 2 J2^ 2nx Pn(x), (2.32) 

n=0 

where the polynomials p n (x) are defined recursively by 

2 d d 2 \ , , 
An + An— + — p n (x) = xp n _i(x) , n = 1, 2, 3, ... , 
ax ax z J 

p (x) = -x - x - 7. (2.33) 

The recurrence relations are spelled out in Appendix iBl along with the proof that the series 
expansion ()2.32j) converges uniformly for x G (— 00, 0). Moreover, we included another 
check that for Nf = the solution given by ()2.32|) reduces to the Bessel function K (pq). 
Indeed, if we substitute y 0jq = 27T 2 K (pq) into ([2.190 . we recover the harmonic function of 
D3-branes, as expected for Nf = 0. This proves that the series (|2.32|) provides the correct 
analytic solution to the differential equation (|2.18|) . 

Upon substituting ()2.32j) in (j2.19j) . we finally obtain the warp factor 

h(r, p) = l + Q m rfg M^j J2 \ n e 2nx p n {x) , x = \og{p/p L ) , (2.34) 

J lT n=0 

where the polynomials p n (x) are given by Eq. ()2.33|) . In Fig. l2.ll we compare the near-core 
solution with the exact solution for y(p). In contrast to the exact solution, the near-core 
solution is only valid for small p and diverges at p^ — e 2n ^ Nf . 
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Figure 2.1: Plot of the near-core solution y(p) = i^o(V 'Xxe 2x ) (x = log p/pi) and the 
exact solution y(p) for Nf = 1, q — 0.1. The near-core solution breaks down around p^. 

2.3 The decoupling limit and regime of validity 

We shall now discuss the regime of validity of the supergravity solution. As mentioned 
earlier, it is in principle possible to find a (numerical) solution to the full D3/D7 configu- 
ration using the complete axion-dilaton of the D7-branes. Our background ()2.1[ 12.21 12.6|) 
approximates this solution in the regime 

0<p<p L . (2.35) 

This regime is bounded from above due to the naked singularity at pi, where the dilaton 
diverges. The other limit is basically p = 0, where the space has a curvature singularity 
due to the presence of the D7-branes. At this point the string coupling goes to zero which 
corresponds to the fact that the dual field theory becomes free in the infrared. 

The decoupling limit is obtained as follows. The massless open string degrees of 
freedom correspond to a Af = 4 super- Yang-Mills multiplet (3-3 strings) coupled to Nf 
bifundamental M = 2 hypermultiplets (3-7 strings) localized at the 3 + 1 dimensional 
intersection. We take a limit in which the 7-7 strings decouple, leaving a purely four- 
dimensional theory. This decoupling is achieved by taking a large N c limit while keeping 
the four- dimensional 't Hooft coupling A = g\ M N c ~ g s N c and Nf fixed. This is the usual 
't Hooft limit for the gauge theory describing the N c D3-branes. The eight-dimensional 
't Hooft coupling A' for the Nf D7-branes is A' = \(2wl s ) A N f / N c which vanishes in the 
a' — > limit. This turns the SU(Nf) gauge group on the D7-branes into a global flavor 
group. 

Moreover, we take the limit of large 't Hooft coupling. Then, the near-horizon ge- 
ometry is obtained, as usual, by dropping the "1" in the warp factor h(r,p). Writing 
g^e* = 1 — uWogp/ pl with v = Nf/N c , one can easily see that in the strict v — > limit, 
the near horizon geometry reduces to AdS$ x S 5 with embedded probe D7-branes pQ. 
Recall that in this limit the warp factor ()2.34|) reduces to the harmonic function of the 
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D3-branes. In other words, the background includes corrections to the quenched ap- 
proximation of first order in Nf/N c . In the field theory these corrections correspond to 
Feynman diagrams with empty quark loops. 

The solution seems to be valid for any number of flavors. Neither the perturbative 
field theory nor the supergravity solution imposes any constraint on the number of flavors. 
Note, however, that a D7-brane solution, which extends to p — > oo, is known to exist only 
for Nf < 12. It is therefore not clear, whether our D3/D7 solution can be extended to a 
full solution in the case of more than 12 flavors. It should be safe to use it for Nf < 12, 
in which case the solution can in principle be numerically continued to infinity. 



2.4 The warp factor of a fractional D3/D7 brane system 

As a bonus for finding the warp factor for the D3/D7 brane system in flat space, we can 
actually construct the warp factor of another system, that of fractional D3/D7 branes on 
the orbifold M 5 ' 1 x C 2 /Z 2 , where the Z 2 action is a reflection parity along the x A , . . . x 7 
coordinates. The D3-branes span the x°, . . . x 3 coordinates and the D7-branes are oriented 
entirely along the orbifold, spanning x°, . . . x 7 . Due to this fact, the D7-branes also source, 
beside the axion-dilaton, a twisted scalar and a twisted 4-form potential. Their charge 
is half that of ordinary D7-branes, for which reason they are called fractional D7 branes. 
The system whose warp factor we will construct next is composed of Nf fractional D7- 
branes and N c fractional D3-branes [HE]- The warp factor has been known so far only in 
terms of the differential equation 

(□ 1 + e-*D 2 )ff + lK 2 „ t T,N c 5(x l )...5(x>) 

where t% is the tension of the fractional D3-brane, denotes the four- dimensional flat 
space Laplacean along the directions {x 4 . . . x 7 } = {r} and \3 2 is the two-dimensional flat 
space Laplacean along the directions {x s ,x 9 } = {p}. Also, the dilaton is 

e^ = l-^lnQ. (2.37) 

We recognize the differential operator acting on the warp factor being the same as the 
one whose Green's function we have computed previously 

G(p, & r; p', <p', f) = £ J d'qe^-^S^'^p^ip^ , (2.38) 

where p — {p, (p}, p< = rnin(p,p'), p > = max(p, p'), and yi, q (p<) is the GP solution, 
well behaved at the origin, and yi, q (p>) is the generalization of the solution discussed in 
App. B. Therefore, the full warp factor of the fractional D3/D7 system is given by 



H(p,r) = 1- J dV J d 2 pG{p,r-J^) 



2Ki b r 3 N c 6\n6\f/) 



\2 



+ Wag.) mi _ Nfg j {2n)HW£)) ** I?) 



(2.39) 
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with the numerical constant C — — 1 fixed by requiring that the discontinuities of Eq. ()2.36j) 
at p = p' match on both sides. It is worth saying that the dual gauge theory of this super- 
gravity background is an d = 4, M = 2 supersymmetric Yang-Mills theory, with SU(N C ) 
gauge group and Nf fundamental hypermultiplets. The orbifold action has projected out 
the adjoint hypermultiplet of the gauge theory dual to an ordinary D3/D7 system, leav- 
ing a gauge theory with a negative beta function for Nf < 2N C . This means that the 
fractional D3/D7 supergravity background with a logarithmically running axion-dilaton 
is dual to the Af = 2 gauge theory at energies above the scale where the gauge coupling 
diverges and where non-perturbative effects become important. 



2.5 The dual M = 2 field theory 

The Af = 2 field theory dual to the near-horizon geometry of the D3/D7 configuration 
in flat space has been studied in [H E]. This theory has U(N C ) gauge group and the 
field content of Af = 4 super Yang-Mills and Nf hypermultiplets in the fundamental 
representation. The global symmetry of the theory is 5*0(4) « SU(2)$ x SU(2)n which 
consists of an SU(2)<& global symmetry rotating the scalars in the adjoint hypermultiplet 
and an SU(2)k M = 2 R-symmetry. In the case of overlapping D3 and D7-branes, the 
fundamentals are massless and the theory has an additional U(l)n R-symmetry. 

We are interested in how far the supergravity background reflects the perturbative 
aspects of the gauge theory. The exact perturbative Af = 2 beta function is given by 

d of 2 o 2 of 2 

(3 N=2 = p— a = -— [2T G - T adj - N f T {und ] = (3 N=4 + —Nf = —N f , (2.40) 

where a = gy M /47r, T G = N c , T ad j = 2T G , T fund = 1 and /?at=4 = 0. In the 't Hooft 
limit, the theory is considered at large iV c with g YM small such that the 't Hooft coupling 
A = g\ M N c remains fixed. In this limit a much more meaningful quantity is the beta 
function of the 't Hooft coupling 

We see that if Nf/N c is kept fixed, then the field theory is neither conformal nor asymp- 
totically free. The theory is however conformal in the strict N c — * oo limit pQ. 

A positive beta function leads to a running gauge coupling gyu which is given by 

1 / 47T \ 

«(Q 2 ) = ~ B ^ with Ai = p 2 exp ——- , (2.42) 

Q 2 the energy scale, p 2 a reference scale and /3 = Nf. A L is the scale of the Landau pole 
at which the gauge coupling becomes divergent. 1 

1 In the N = 2 theory of the orbifolded D3/D7 intersection considered in Sec. 12.41 adjoint hypermul- 
tiplets are projected out, and the beta function (3 X ~ —(2 — Nf/N c ) is negative, if Nf < 2N C . 
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This can be compared with the action for a probe D3-brane placed in the D3/D7 
supergravity background (|2.1[ \'2.'2\ I2.6jl . 

S D3 = ~T m j d 4 ae-*y/-det(g ab + F ab )+ T D3 J C A + C°F A T 



D3 



d 4 a(2iia 



'\2 



-l e -*F ab F ah + X F ab F ab 



+ ..., (2.43) 



where T ab = (2iia')F ab . Here we expanded to second order in the field strength and kept 
only terms relevant for the field theory. The axion C° = x an d the dilaton are given by 
the complex structure (jA.5|) . They are 



x =%, e-* = ^log4. (2.44) 

27T 47T p 2 



The D3-brane action relates the dilaton and the axion to the gauge coupling g\ M = 
47re* and theta angle 9ym — ^X- Upon identifying also a(Q 2 ) = g Y MiQ 2 ) /47r, a(p 2 ) = 
g YM (p 2 )/4:iT = g s , Q 2 = p 2 /(27ra' 2 ), A L = p L /(2na' 2 ), the running of the gauge cou- 
pling (J2.42j) follows from the logarithmic behavior of the dilaton, and the Landau pole 
{qym oo) is related to the dilaton divergence at 



p L = e 9sN i =e XN f . (2.45) 

This shows that both sides of the duality show the same pathology: The perturbative field 
theory becomes strongly coupled at the Landau pole A L , while the supergravity solution 
breaks down at some distance pi 2 

Furthermore, the U(l)n chiral anomaly in the field theory is reflected by a nontrivial 
axion profile in the supergravity background. The instanton term in ()2.43j) correctly 
reproduces the Yang-Mills theta angle 9ym = — Nfip. This corresponds to the 
breaking of the classical U(l)n symmetry to Z 2 N f at the quantum level. 

This shows that the solution (j2.1j 12 . 2| 12. 6|) of the D3/D7 system perfectly reflects the 
perturbative aspects of the dual Af = 2 field theory. 



2.6 Logarithmic tadpoles and one-loop vacuum amplitudes 

Since D7-branes are co dimension- two branes, we expect uncanceled tadpoles in the string 
background. Tadpole divergences usually correspond to gauge anomalies and indicate an 
inconsistency in the theory. However, as it was found first in logarithmic tadpoles 
do not correspond to gauge anomalies, but reflect the fact that the dual gauge theory is 
not conformal invariant. In fact, as we will show now following arguments given in |34j . 
such tadpoles provide the correct one-loop running of the gauge coupling. 

This can be seen by considering the one-loop annulus diagram of an open string stretch- 
ing between the stack of D7-branes and a D3-brane dressed with an external background 

2 In the conformal limit Nf/N c — > 0, the location is formally shifted to infinity and the near-horizon 
geometry of the solution reduces to AdS 5 x S 5 with embedded probe D7-branes. 
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SU(N C ) gauge field. In the field theory limit this vacuum amplitude reduces to the one- 
loop diagrams determining the running of the gauge coupling [HI]. When expanded to 
quadratic order in the gauge field, the amplitude reads: 3 



Nt 



J 37 



(4tt) 



d 4 x 



4 K r nu r ir tiv r 



dr 
— ( 

l/(a'A 2 ) r 



(2.46) 



The integral over r is logarithmically divergent for small r and has been regularized by 
an ultraviolet cut-off A. The parameter p is the distance between the dressed D3-brane 
and the D7-branes. Threshold corrections to the amplitude ()2.46|) . coming from massive 
string states, are absent in the a' — > limit. 

From the annulus amplitude, one deduces the gauge coupling 



9ym 



(An) 



dr 



N f 



l/(a'A 2 ) 



T 



(47T) 



1 A2 

log g^' 



[2A7) 



where we defined the energy scale Q 2 = p 2 /(2ira' 2 ). The previous relation agrees with 
Eq. ()2.42|) and the dilaton behavior described by ([2.44)1 . This shows that tadpoles provide 
the correct one-loop running of the gauge coupling. 

Furthermore, from the instanton term in (J2.46j) . we may also extract the 6ym angle, 







YM 



Nfip. 



(2.48) 



where p is the phase of the complex coordinate z = pe tip transverse to the D7-branes. As 
mentioned above, in the supergravity background the chiral anomaly is encoded in the 
axion x- 

Under open/closed string duality, one can consider the above one-loop amplitude also 
as a tree-level closed string amplitude which encodes information about the supergravity 
solution. The absence of threshold corrections, i.e. massive open string state contributions 
to the gauge coupling, guarantees that open massless string states are precisely mapped 
into massless closed string states. The supergravity solution constructed in this paper is 
therefore sufficient to describe the perturbative field theory. Even though the dual string 
background has logarithmic divergences, it is still consistent as far as the description of 
the embedded non-conformal field theory is concerned. 

Finally, it should be clear that in the conformal limit Nf/N c —>■ 0, considered by Karch 
and Katz pQ , the gauge coupling does not run, and by the above arguments, tadpoles are 
absent. Certainly, D7-branes which do not backreact on the geometry do not emit flux 
and there is no net charge. 



3 This one-loop open string amplitude is obtained by summing up the one-loop amplitudes of the 
twisted and untwisted sector of the orbifolded D3/D7 configuration which have been computed in |34| : 
it is identical to 2Z e: 3i = {Z e -^j + Zh-.zi) + (-Ze;37 — Zh-,3T) given by Eq. (75) therein. Note also that the 
amplitude is different from the one corresponding to a string stretching between the stack of D3-branes 
and the stack of D7-branes. The latter vanishes as expected for a BPS Dp-D(p+4) brane configuration. 
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Figure 3.1: Example of a string profile z(R). 

3 Regge trajectories in the D3/D7 theory 

As an application of our supergravity solution, we now discuss the general structure of 
Regge trajectories in the M = 2 gauge theory of the D3/D7 system by means of a semi- 
classical string computation. This will provide interesting information about the behavior 
of the Regge trajectories in dependence of the number of flavors. We shall compare our 
results with a similar analysis performed in the probe approximation in [2j. 

3.1 Spinning strings in the D3/D7 background 

Following [SHI 121 E]; we consider an open string rotating in the near-horizon limit of 
the D3/D7 background. We parameterize the four- dimensional spacetime part of our 
background (EH O EH as 

dx^dXf, = -dt 2 + dR 2 + R 2 d</ + dz 2 , (3.1) 

where R and <p are the coordinates of the plane of rotation. The string has length 2Rq 
and stretches from — Rq to +Rq along the R direction. The end points of the string 
are attached to an additional probe D7-brane located a distance pR away from the stack 
of D7-branes. Recall that we have to choose pu in the regime p# p^, i.e. far 
away from the singularity at the location of the stack of D7-branes and from the dilaton 
divergence at pi . An example of a spinning string is shown in Fig. 13.11 

In the field theory the set-up corresponds to Nf massless flavors plus an additional 
massive flavor whose mass is proportional to p#. In this theory we consider a meson, which 
consists of massive flavors, in the presence of Nf massless flavors. We assume a large spin 
for the meson which allows a classical treatment of the dual spinning string. Recall 
that meson operators with large spin have small anomalous dimensions and quantum 
corrections are negligible [3Hj . 

An appropriate ansatz for a string rotating with constant angular velocity uj is 

t = r, ip = uiT, R = R(a), r = r(a), p = p(a) , (3.2) 
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with world sheet coordinates a and r. With this ansatz the classical Nambu-Goto action 
takes the form 



C = -T 8 \H\- ujm^ih- 1 !}' 2 + r' 2 + (1 - ^ log p)p' 2 ) . (3.3) 
V 2ir 

It is convenient to use the rescaled coordinates 

R = ujR, f=-, p = -. (3.4) 

UJ UJ 

In these coordinates, the energy and the angular momentum of the spinning string are 
given by 

E= I do(uj^--C \ = I dajyJh-iR'ia) 2 + f> (a) 2 + e^p'(a) 2 , (3.5) 



dc r n 



J= jda—= da—sJh^R'^Y + f^af + e^p^a) 2 , (3.6) 



with £ = Vl — R 2 , T s — 1 and 



^^ )= (^f' (3,7) 



In the gauge R = a, we find the following equations of motion for r(R) and p(R): 

d (£\ A 



(9j=,r ) = -^dfh 1 , (3. 



The warp factor h(f, p) satisfies dfhr 1 = for f = which shows that f = is a solution 
of the equation of motion ([3.8)1 . Writing Eq. (|3.9|) in coordinates 5 = 1/p, we find 

This is a nonlinear differential equation of second order which requires two boundary 
conditions. We impose the usual open string boundary conditions 



™6R 
dR' 



= (3.11) 

(7=0, 7T 



CT=0,7T 

for a string ending on a probe D7-brane at z = zr = const. Due to the Neumann boundary 
condition in the R direction and the Dirichlet boundary condition in the z direction, 
5R\o-=o,n is arbitrary, whereas 8z\ c= q i% = 0. The remaining condition dC/dR'\ a= o,n = 
is satisfied, if i?'| cr=0i7r = 0. Using the gauge z = a, we see that this corresponds to 
dz/dR\ R=±Ro — > oo. This means that the string ends orthogonally on zr. 

We expect the solutions to be symmetric around R = 0, where they have their only 
maximum. We thus impose the additional boundary condition z'(0) = 0. Taking into 
account the orthogonal ending of the string on a constant value of z, we set z(±R ) = 
zr = const, or, equivalently, 5(0) = z = const. 
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Figure 3.2: Chew-Frautschi plot for Nf = 0, 1, 2, 3, 5, 10 additional massless flavors. The 
straight line represents the Nf = trajectory for small spin values. All graphs approach 
the horizontal line E = 2m. 

3.2 Dependence of Regge trajectories on the number of flavors 

The Regge trajectories E( J; Nf) can now be obtained as follows. We first solve the 
equation of motion (|3.1Uj) for the string profile z(R). We then substitute the profile into 
the expressions (|3.5|) and (|3.6|) for the energy E(u) and the angular momentum J(uj) 
of the spinning string. This yields a Regge trajectory as a curve in the y/J, E plane 
parameterized by the angular velocity uj. 

We begin by fixing the mass m of the single flavor corresponding to a 3-7 string 
stretching from the D3-branes at p = (z — > oo) to the probe D7-brane at zr = uj/pr. 



where e <C 1 is a (IR) cut-off which is needed to regularize the singularity at p = 0. We 
set m = 1 and solve this equation numerically for pn(m, Nf). In this way we obtain the 
location of the probe D7-brane ZR{m, Nf) = uj / pR{m, Nf). The set of Regge trajectories 
in Fig. 13.21 is obtained by varying u for a given number of massless flavors Nf, then 
changing Nf while all the time keeping the mass m fixed. 

We now determine the string profile zr for given Nf and u. To this end, we integrate 
the equations of motion (j3.10j) from — R Q to +Ro- In the shooting technique, we set 
z(0) = z Q = const, z'(Q) = such that z(±R ) = z R . This yields the string length 
Rq = Rq/oj as the location at which z'(Rq) — > oo. A typical profile has been shown in 



4 In general, one obtains a series of solutions labeled by the number of nodes of the string We 
restrict to string solutions with no nodes which are believed to be the most stable ones. 



m = 




(3.12) 



Fig. mm 4 
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Figure 3.3: String length Ro(J) for different Nf. 

Substituting the profiles r(R) = and p(R)~ l = z(R) into Eqns. (j3.5|) and (j3.6J) . we 
determine the energy E(uj) and the angular momentum J(uj) in dependence of the angular 
velocity uj. Repeating the procedure for different u>, we obtain a parametric plot for the 
Regge trajectory E(\/j; Nf) associated with a certain number of flavors Nf. 

In the computation we keep g s ~ A/iV c small, but fixed such that the background 
depends on the small but finite parameter v = Nf/N c . 

Our results are presented in Fig. 13.21 which shows the Regge trajectories for different 
numbers of massless flavors Nf. 5 The energy E and the spin J of the mesons are given 
in units of the quark mass m and the square root of the 't Hooft coupling A, respectively. 
As a first check of our numerics, we note that the graph for Nf = coincides with that 
found in the probe limit 0. Disregarding the behavior at small spin values, the graph 
for Nf = runs above the graphs for Nf ^ 0. 

The semiclassical approximation is valid as long as the quantum effects are suppressed 
by T s /S c is <C 1 where S c i s = E/u — J is the value of the action evaluated on the classical 
solution. One can check that while the linear behavior of the Regge trajectories lies well 
outside the region of validity of the semiclassical approximation, 6 the region where the 
quantum effects are suppressed is still sensitive to the number of flavors. 

Let us first consider the general behavior of the Regge trajectories. Quite generically, 
the trajectories are linear in the small spin region and asymptote to the rest energy 
E = 2m in the large spin limit. This can be understood from the behavior of the string 
length as a function of the spin. Fig. 13.31 shows that the data fits a linear relation between 

5 For each graph there exists a spin value J at which the location zr of the probe D7-brane becomes 
of the order of zl- Here the solution 1)2.11 12. 2| 12.6(1 breaks down. The dashed part of the graphs is an 
interpolation to higher spin values. 

6 Thc spinning string characterized by small E, J is not macroscopic and probes a small region of 
the warped geometry which, being small, appears almost flat giving rise to linear Regge trajectories. 
Quantum effects however will largely correct this apparent confining behavior. 
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Figure 3.4: String tension in dependence of the number of flavors Nf for J = OA 2 (upper 
graph) and J = 0.6 2 (lower graph). 

the string length Rq and the spin J. Note that this plot contains data from graphs for 
different Nf. Since the linear relation is the same for all Nf, we find that the string length 
Ro does not depend on the number of massless flavors. 

At small spin values the length of the string is much smaller than the scale of the space, 
-Ro Qd3, and the string is effectively rotating in flat space leading to a linear Regge 
behavior. At large spin the string is larger than the size of the space, Rq 3> Qd3- Here 
the string rotates very slowly and the energy is that of particles moving in a Coulomb 
potential 0. Moreover, since the string length grows linearly with spin, also the binding 
energy of the quark- ant iquark pair vanishes at large spin values. The Regge trajectories 
thus asymptote to the rest energy E = 2m in the limit uo — > 0. 

Let us now compare trajectories with different numbers of flavors. We observe that for 
a meson with fixed (large) spin, the energy is the lower the larger the number of flavors in 
the theory, i.e. E\ > E2 for Nj < Nj and J\ = J^. In other words, the graph for smaller 

Nf is above the graph for higher Nf in the \/~J, E plane. 

To understand this behavior, it is illuminating to consider the string tension for a 
fixed spin value as a function of Nf. The square root of the tension is proportional to the 
derivative of E{\fj). Fig. 13.41 shows that the string tension is a monotonically decreasing 
function of Nf. Since other quantities like the string length are unaffected by an increasing 
number of flavors, the energy of a meson is smaller in a theory with more flavors. This 
effect vanishes however at very large values of the spin. 

The behavior of the string tension can also be understood in the dual field theory. 
The rotating string can be considered as a spin chain of gluons with quarks at the end 
of the string. In between the quark and the antiquark of the meson there are virtual 
quark- ant iquark pairs which interact with the gluon chain. The polarization of these 
virtual quark-antiquark pairs effectively diminishes the color force between the quarks 
of the meson. This screening effect should have an influence on the string tension. We 
expect that a large number of flavors in a field theory leads to a small string tension. At 
least heuristically, this explains the fact that a meson with a certain spin has less energy 
than the same meson in a theory with a larger number of flavors. 
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Even though the background we are discussing is not confining, a similar effect is 
observed in lattice QCD with dynamical light quarks, where the interaction potential 
between two heavy quarks 

V(r)=( + ar) (3.13) 

\ r J pr 

was a better fit to the data, rather than the unscreened potential jlU] . For non- vanishing 
inverse screening length p, the effective flux tube tension between the heavy quarks is 
diminished due to the screening of the color charge, with p increasing with the number 
of light quarks. 



4 Conclusions 

To summarize our results, we have completed the type IIB supergravity solution describing 
the fully localized D3/D7 brane configuration by providing an analytic expression for the 
warp factor. The background is given by the metric 

ds 2 10 = /T 1/2 (r,p) dx\ + /i 1/2 (r,p) (dr 2 + r 2 dn 2 3 + e^ p \dp 2 + p 2 d<p)) (4.1) 

for which we assumed a logarithmically running axion-dilaton: 

AT, AT, n 2tt 

X(f) = > e ~ m = e = -IT log ' PL = e93Nf ■ ^ 
An In Pl 

Upon taking the decoupling limit, this is the supergravity dual of M = 4 SU(N C ) super 
Yang-Mills gauge theory with Nf M = 2 fundamental hypermultiplets. 
The warp factor h(r, p) is given in terms of a convergent series: 

h(r, p) = l + Q m f°° dq {qr)2Jl } qr) V \ n e 2nx Vn {x) , x = log(p/p L ) , (4.3) 

where A = -^-PlQ 2 , and the polynomials p n (x) (n = 0, 1, 2, ...) are defined recursively by 
d d 2 \ 

An 2 + An— + — p n = xp n _i , p (x) = -x - \og(p L q/2) - 7 . (4.4) 
dx dx z J 

In the immediate vicinity of the D7-branes, x — > —00 (p — > 0), the warp factor can be 
approximated by the simple expression 

Kr,p) = i+ (ra+ ^ M)a - (4.5) 

We have related the pathology of the supergravity background to the Landau pole 
of the gauge theory which has a positive beta function. We have also shown that the 
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existence of a logarithmic tadpole in the one-loop open string amplitude between the D3 
and D7-branes is related to the same one-loop running of the gauge coupling. 

In the second part of this paper, we have considered classical spinning string configu- 
rations in the backreacted D3/D7 geometry and analyzed the dependence of the associ- 
ated meson Regge trajectories on the number of flavors. From the Chew-Frautschi plots 
(E, vJ) we have found that for the same spin, the energy of a meson is smaller in a gauge 
theory with more flavors. We interpreted this in terms of a screening effect of the color 
charge due to virtual qq pairs. Since the background is non-confining, in the IR, where the 
theory becomes conformal, the tension of the quark-antiquark flux tube vanishes at large 
spin J, and the meson configuration energy approaches the rest mass of its constituents. 

We have seen that taking into account the backreaction of the flavor branes opens 
up the possibility of going beyond the quenched approximation via the AdS/CFT corre- 
spondence. An immediate goal would be the study of confining backgrounds, including 
the backreaction of the flavor branes. Also, it is well-known that SU(N) gauge theories 
undergo a chiral phase transition [22] as the number of flavors is varied, even at zero 
temperature. It would be interesting if these phase transitions could be studied from the 
perspective of the dual supergravity background. 
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A Brief review of the D7-branes geometry 

For completeness, we review here the geometry of a stack of D7-branes which has first 
been studied in jjl] in the context of cosmic string solutions. The standard metric ansatz 
for the type IIB supergravity solution corresponding to a stack of Nf coincident D7-branes 
is 

dsl = dsl + e t ^dzdz, (A.l) 

where ds^ is an eight-dimensional flat metric and z is a complex coordinate parameterizing 
the space transverse to the D7-branes. D7-branes couple to the axion-dilaton field r = 

X + ie~^. 

The complex structure r and the function ^ have to satisfy the equations of motion [JT] 

— 2c)t8t 

ddr + — — = N f S(z) , (A.2) 
r — r 

dM = dd log Im (r) . (A.3) 
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The first equation requires r to be holomorphic or anti-holomorphic. In order to have a 
finite string coupling, Im (r) must be nonzero for all values of \z\. The basic solution for r 
can be written in terms of the modular invariant j-function A stack of Nf D7-branes 
corresponds to a pole of order Nf in the j-plane which can locally be parameterized by 

J(r)= 1 (A.4) 

{z/pL) Nt 

For small \z\, j(r) is large and can be approximated by exp(— 2nir). Substituting j ~ 
exp(— 2Trir) in Eq. (|A.4|) and solving for r, we get the weak-coupling approximation 7 

.N f . fz 



^—^{vj- (A - 5) 



The integration constant pi determines the string coupling. 
The second equation of motion, Eq. (jA.3J) is solved by 

e^ z ^ = Qtt Im (r) , (A.6) 

where Q is an arbitrary holomorphic function. This function is fixed by requiring modular 
invariance and regularity of the metric at the location of the D7-branes. The function e* 
is thus given by 

e 9 ^ = t 2 (z)\ V (t)\ a \z\- n ^ (A.7) 

with 7] the Dedekind eta function and r 2 = Im (r). In the weak coupling region, \z\ <C pl, 
t 2 is large and the function e* is well-approximated by 

e* {z ~ z) w r 2 . (A.8) 



Asymptotically, for \z\ 3> Pl, the complex structure r as given by Eq. (jA.4|) approaches 
the constant r = j _1 (0) and e* <z ' 2 - ) ~ l^l"^/ 6 . This leads to an asymptotically flat metric 
with a deficit angle of 2irNf/12 as can be seen from the transverse part of the D7 metric, 

dsl = p- N ^(dp 2 + p 2 V) = — l -WT 2 W 2 + P' 2 (l - V) , (A.9) 

^ 12 > 

where p' = p l ~ N fl 12 . This shows that the transverse space is asymptotically conical for 
Nf < 12 and asymptotically cylindrical for Nf = 12. For Nf > 12 the transverse space is 
compact and in general singular apart from the exceptional case Nf = 24. 

7 This is the complex structure which one would naively expect from the general behavior of the dilaton 
of a Dp-brane (p < 7), = (1 + A/ p 7-P)( P -3)/4_ 
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B Solving the equation (12.181) 

The second order differential equation (|2.18|) admits two independent solutions. One of 
them, which is well behaved near x = 0, was found by Gesztesy and Pittner 



y(x) = J2^ 2nX Pn(x), (B.l) 



n=0 

where the polynomials p n (x) are defined by the recursion relation 

An 2 + An^ + -^^j p n (x) = xp n -i(x) , p (x) = 1 . (B.2) 

The series converges uniformly on the negative real line, x < 0. 

Here we proceed to derive the other independent solution which is needed in the 
construction of the warp factor of the backreacted D3/D7 geometry. 

The recurrence relation 

An 2 + An^- + ^)pn{ x ) = xp n -i(x) (B.3) 

can be decomposed into a set of algebraic equations for the expansion coefficients of the 
polynomials p n (x). With the initial condition 

Pq{x) = ax + b , (B.4) 

where a, b are constants, we find that 

n+l 

Pn(x) = ^4 n) x fc , (B.5) 

k=0 

and where the coefficients c^. satisfy the following system of equations 

An 2 c^ + An(n + = , 

An 2 c { :l k + An(n - k + l)c^ fc+1 + (n - k + 2)(n - k + l)c^l k+2 = k > 1 

4n 2 4 n) +4ncS n) + 24 n) = 0. (B.6) 

According to our initial condition, we have = b, cf = a. One can easily solve for the 
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leading expansion coefficients (for x — > — oo) 



» 

--n+l 



» 
-n-1 



» 



4™(n!) 2 ' n 4™(n!) 2 4™(n!) 2 

_J_( 4 „ 2c ( r ) +n(n+1)cSi+ i 



(n + ^(n+ 1) +7) 



4(^-l) 2 c!r i 1) +n(n-l)ci"- 1 ) + 



4(n-2) 2 



4(n - 2) 2 cir 2 2) + (n - l)(n - 2)^ + . . ^ + 1 • 2c« 

ta _ fl |^ (n+ l) +7) + ^ " 



4"(n!) s 



(B.7) 



where ip(n) = ^"d}^ ls ^ e polygamma function, and 7 = is the Euler-Mascheroni 

constant. 

Let us show that in the limit when the number of flavors is set to zero (i.e. there 
are no D7-branes), we recover the modified Bessel function K (qp) which is needed to 
reconstruct the warp factor of the geometry obtained by taking the near-horizon limit of 
a set of overlapping D3-branes (i.e. the Green's function of the six-dimensional flat space 
Laplace equation). 

We begin by reintroducing the initial set of variables 



2tt 



N 



f 



(B.8) 



which yield 

y(x) = 5>e 2 *)>n(*) 



+ ( b - a ( n + 7 + i>(n + 1) ) ) (In p - ^-) r 



+ [bn-a(n + 1) 
+ 



(2±^I +7 + Wn + 2 )))(h.p-^) 



In order to obtain 



^ (gp) 2 "(^(n + l)-log(gp/2)) 

^o(gp) = ^ 4^ , 



n=0 



we must choose 



27ra , . . . 

& = -rr- + cry + aln(g/2) . 



(B.9) 



(B.10) 



(B.11) 
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Notice that with this choice we find that the only terms in the x(p) expansion of y(x) 
that contribute in the limit Nf — > are at most of order x n , and at this order, the only 
contribution comes from the 6-dependent term. The anN^x n term is now canceled upon 
taking Nf — > by the bnNJx n ~ 1 term. This completes the proof that y(x(p)) reduces 
indeed to Ko(qp) when setting the number of flavors to zero. Moreover, the requirement 
that in this limit we obtain the modified Bessel function K Q (qp), allowed us to enforce the 
relation (jB.ll|) between the ab initio free parameters a, b. The remaining dependence on 
a manifests as an overall proportionality coefficient which we fix by the same requirement 
that in the limit Nf — > the solution y(x(p)) becomes K (qp). From now on, we set 
a = -1 and (lB~TT|) . 

The absolute value of the coefficients is bounded from above by 

l4 n) l = (\a\Mn + 1) + 7 + l)^}-*^ 1 "* + |&|CT*2»-*) • (B.12) 
Consequently, the polynomials p n (x) obey the inequality 

\Pn(x)\ = JL7|a|(V>(n + 1) + 7 + 1)(2 + M)" +1 + |6|(2 + . (B.13) 

Therefore it can be shown that the series y(x) = X n e 2nx p n (x) is convergent given that 
each term in the series is smaller than the corresponding term of a convergent series 

E \M ne2nx ^.\\ a \(i ) ( n + 1) + 7 + 1)(2 + + |6|(2 + The convergence of 

the latter series is proven by showing that the ratio of two consecutive terms a„/a n _i = 
Ae x (2 + x)(|a|(2+|x|)(^H + 7+l) + |&|)/((a(2 + |x|)(^(n-l)+7 + l) + |6|)4n) approaches 
zero as n — > oo. 
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